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Abstract 

A Pick function of d variables is a holomorphic map from U d to II, 
where II is the upper halfplane. Some Pick functions of one variable 
have an asymptotic expansion at infinity, a power series Y^=\ Pn z ~ n 
with real numbers p n that gives an asymptotic expansion on non- 
tangential approach regions to infinity. H. Hamburger in 1921 char- 
acterized which sequences {p n } can occur. We give an extension of 
Hamburger's results to Pick functions of two variables. 

1 Introduction 

A Pick function of one variable is a holomorphic map from the upper half- 
plane, which we shall denote by IT, into II. A Pick function of two variables 
is a holomorphic map from II 2 to II. The purpose of this note is to extend 
to two variables certain well-known results about the asymptotic analysis of 
Pick functions in one variable. 
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1.1 One variable results 



In 1922, R. Nevanlinna showed that a Pick class function of one variable that 
decays at infinity is the Cauchy transform of a finite measure on R. 

Theorem 1.1. [13] If F : U — > II is analytic and satisfies 

limsup \yF(iy)\ < oo, (1.2) 

then there exists a unique finite positive Borel measure p onM so that 

F W = / (1.3) 

nt 

We shall say that a set S in II approaches oo non-tangentially, S — > oo, 
if oo is in the closure, and there is a constant c such that \z\ < c Im(z) for 
all z G S 1 . If .F has a representation as in (1.3), then 

F{z) = P - + o{l/\z\) 

as z ^> oo, where p = — If /i has more moments, then there is a higher 
order asymptotic expansion at oo. H. Hamburger proved the following two 
theorems [8, 9]. For a proof of Theorem 1.4 as stated, see [17, Thm. 2.2] or 
[5, Thm 3.2.1]. 

Theorem 1.4. Let real constants p±, . . . , P2N-1 be given. There exists a Pick 
function F satisfying 

F{z) = ^ + ^ + --- + ^ + o(N" (2 "" 1) ) (1-5) 

as 2 4 00 i/ and only if there is a measure p on R whose first 2(N — 1) 
moments are finite and satisfy 

J t k dp{t) = - Pk+1 , < k < 2(N - 1). (1.6) 

Moreover, in this case F has a representation as in (1.3) for some measure 
p satisfying (1.6). 
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Hamburger gave an alternate equivalent condition. There is also a proof 
in [17, Thm. 1.2]; and see [12, Thm. 3.3] for an alternative formulation (but 
without a proof). 

Theorem 1.7. Let pi, . . . ,P2N-i be given real numbers. There exists a Pick 
function F satisfying (1.5) if and only if the N-by-N Hankel matrix 

'Pi P2 ■ ■ ■ Pn \ 

Pi P3 PN+1 

rl — — 

\ Pn Pn+i ■ ■ ■ P2N-1 J 

is positive semi- definite and has the property that whenever {c\, C2, . . . , cjv-i, 0)* 
is in the kernel of H , then (0, ci, . . . , cn-i) 1 is also in the kernel. 

In 1881, L. Kronecker proved the following theorem [11] (see [14, Thm. 
1.3.1] for a modern treatment). 

Theorem 1.8. The infinite Hankel form 

( Pi P2 P3 
P2 P3 P4 
P3 P4 P5 

\ ; : ; 

is finite rank if and only if 

00 

*w = 

n=l 

is a rational function. 

1.2 Two variable results 

A two variable version of Theorem 1.1 was proved in [4]; see also Theorem 8.4 
below. Before stating it, let us introduce some notation. If Y is an operator 
on a Hilbert space, and z — (z\, z%) is a point in C 2 , we shall use zy to denote 
the operator 

z Y = Zl Y + z 2 (I-Y). 
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Theorem 1.9. [4] Let h : U 2 — > II be a Pick function of two variables. Then 



limsup \sF(is, is)\ < oo 



if and only if there is a Hilbert space H, a self-adjoint densely defined operator 
A on H, a positive contraction Y on H, and a vector a in %, such that 



We shall say that h has a type I Nevanlinna representation if it has a 
representation as in (1.10). 

In one variable, the Poisson integral of any finite positive measure on IR 
is the real part of a Pick function that decays like (1.2), so the study of 
asymptotic expansions (1.5) and solutions to the moment problem (1.6) for 
arbitrary measures are tightly bound. In two variables, their study diverges. 
The infinite Hamburger moment problem in several variables is studied in [16] 
and [18]; for an algorithm for solving the problem in two variables, see [19]. 
For the truncated problem, see for example the memoir [7] and subsequent 
papers. Our objective is to study the two variable analogue of (1.5). 

If one restricts z to the diagonal {z± = z 2 }, then (1.10) becomes (1.3), 
where fi is the scalar spectral measure of A for the vector a. Saying that an 
even moment j2k exists in this case is the assertion that t h ~ l is in the domain 
of A. We shall generalize this idea to two variables. 

We shall let m and n denote ordered pairs of nonnegative integers. We 
set ei = (1, 0) and e 2 = (0, 1). If n = (ni, n 2 ), we set \n\— ni + n 2 , and for a 
pair z = (zi,Z2) we follow the usual convention of letting z n = z^z^ 2 . For 
N a positive integer we set In = {n | 1 <|n|< iV}. 

We now define an object that we shall call a finite Hankel vector moment 
sequence, or for short, a finite HVMS. For simplicity, we take N > 2; see 
(2.1) for general N. 

Definition 1.11. For a fixed positive integer N > 2, a finite Hankel vector 
moment sequence is a 3-tuple, ({a n } nE ; w , Y, A) where: {oc n } neI is a se- 
quence of vectors in some Hilbert space H; Y is a positive contraction acting 
on Ti, satisfying for each I = 1, . . . , N 



A is a partially defined symmetric operator on H, with the property that 



h{z) =<{A-Zy) 1 a,a>, z G II 2 . 



(1.10) 



Ya w) = = (1 - Y)a {tfi) = 0; 



(1.12) 



{«„ 1 1 < |n| < N - 1} C Dom(A); 



(1.13) 



4 



for each n G In-i> 

Aa n = Ya n+ei + (1 - Y)a n+e2 . (1.14) 
Here is the main result of this paper. 
Theorem 1.15. A Pick function h of two variables satisfies 

h(z) = £ £ + oi\\z\\-W-V) (1.16) 

as z — > oo, /or some real numbers p n , if and only if it has a representation as 
in (1.10) and there is a finite HVMS ({a„} ng/iv , Y, A) with a = a(i,o) +«(o,i) ■ 
Moreover, pk is given by the formula: 

pk = ~ ^2{(a n ,Aa rn ) : m 1 + n 1 = fci, m 2 +n 2 = k 2 , m 1 +m 2 = [\k\/2\ }. 
When k — 1, one interprets the right-hand side of the inner product as a 

Tit 

(so P(i,o) = — (tt(i,o)) a ) an d P(o,i) — —( a (o,i)i a ) )• By z — > oo we mean that 
||z|| —7- oo while z stays in an approach region 

{z G II 2 : ||2;|| < c min-flm^!, Im2 2 }} 

for some c. The notation |_7V^T/2J stands for the greatest integer less than or 
equal to M/2. 

The forward implication of (1.15) is Theorem 4.2; the converse is The- 
orem 3.10. To relate Theorem 1.15 to Theorems 1.4 and 1.7, think in one 
variable of ct n as t n ~ x in L 2 (p), and A as multiplication by t on L 2 (p). Then 
Pk is given by a single term, — (t^/ 2 !" 1 , tL fc / 2 J^ 

Theorem 1.7 also has a two variable analogue, which we give in Theo- 
rem 5.7. This justifies our nomenclature of Hankel vector moment sequence. 
The last condition in Theorem 5.7 is an analogue of the last condition in 
Theorem 1.7; for an explanation of it, see Section 5. 

Theorem 5.7 Let a = (a 1 , a 2 ) be a pair of matrices on 1^. Then there 
is a finite HVMS ({a n } n6/jv , Y, A) such that 

a mn = (Y a m a m) 

a mn = ((1 Y)(y n , Ot m ) 
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if and only if the following four conditions obtain: 

a 1 and a 2 are positive semi- definite. 

t n+e 2 whenever m, n E ijv-i- 
a (o,o,(o,0 = a («,o),(z,o) = f or I = I, ■ ■ ■ , N. 
supp(/) E 7jv-i and (a 1 + a 2 )/ = =>- (a 1 S 1 + a 2 5 2 )/ = 0. 

In Section 6, we discuss infinite sequences. One multi- variable gener- 
alization of Kronecker's Theorem 1.8 was proved by S. C. Power [15]. In 
Theorem 6.6, we prove another. 

Theorem 6.6: Let h have non-tangential asymptotic expansions of all 
orders at infinity. Then there is an infinite HVMS ({a n },Y, A) with a = 
a (i,o) + Q!(o,i), an d h( z ) =< 0*4 — zy)~ 1 a,a >. The sequence can be chosen 
with rank(a n , a m ) < oo if and only if h is a rational function. 

In Section 7, we give an example of a construction of functions in the Pick 
class that have asymptotic expansions. In Section 8, we give some technical 
results on models. 

2 Finite Hankel Vector Moment Sequences 

Definition 2.1. For a fixed positive integer N, a finite Hankel vector mo- 
ment sequence is a 3-tuple, ({a n } n <- lN ,Y, A) where: {a n } nelN is a sequence 
of vectors in some Hilbert space %; Y is a positive contraction acting on %, 
satisfying for each I — 1, ... ,N 

Ya m = = (1 - Y)a m = 0; (2.2) 

A is a partially defined symmetric operator on % with the properties that, if 
N>2 

{a n \ 1 < \n\ < N- 1} C Dom(A); (2.3) 
and, for each n E ijv-i, 

Aa n = Ya n+ei + (1 - Y)a n+e2 . (2.4) 

When N = 1 conditions (2.3) and (2.4) are dropped. 
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Every symmetric operator has a self-adjoint extension on a possibly larger 
Hilbert space; so there is no loss in generality in assuming A is self-adjoint. 

If ({a n } n€ i N , y, A) is a finite HVMS, we frequently shall abuse the nota- 
tion somewhat and refer to the entire tuple by simply {a n }. If {a n } is an 
HVMS as above, we refer to N as the size of {a n }, Y as the Hankel weight 
of {o n }, A as the Hankel shift of {o n }, and finally the vectors, a n are called 
the vector moments of {a„}. 

Our first proposition gives a simple yet fundamental property of HVMS's. 
If z G C 2 and Y is a positive contraction on a Hilbert space "H, we defined 
zy = z{Y + z 2 {l — Y). As Y is a positive contraction, the spectral theorem 
implies that Zy IS cL well defined analytic operator valued function on the 
set {zeC 2 \z 2 ^ 0, Zi/z 2 i (-00, 0]}. If {«„} is an HVMS with shift A and 
weight Y, and / is a positive integer we shall adopt the notation, 

R t (z) =z Y \Az Y 1 ) 1 - 1 . 

Note that if z G {z G C 2 | 22 7^ 0, £1/2:2 ^ (—00, 0]}, then the domain of Ri(z) 
is all of % if I = 1 and for / > 2 is inductively defined by 

Dom( (#;(»)) = {a G "H I (z^A/z^a G Dom(A) z = 0, 2}. 

Note also that 

Proposition 2.5. Lei {a n } be an HVMS of size N and let 

ot = 01(1,0) + a(o,i) • (2.6) 

Ifl<l<N, then 

a G Dom^z)) (2.7) 

and 

i? z (z)a = ^ — a n (2.8) 

|n|=i 

/or a// z in {z G C 2 | z 2 7^ 0, ^1/^2 ^ ( — 00 > 0]}. 

Proof. We induct on V. If iV = 1 and 1 = 1, then trivially 2.7 holds. Also, 
by 2.2, Vo (0 ,i) = = (1 - Y)a {1>0) = 0. Hence, 

R\(z)a = Zy l a 
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Zy 0(1,0) + Zy Q!(0,1) 

1 1 

— 0(1,0) + —O(i,o) 



z 

\n\=l 

Now assume that the proposition holds for HVMS's of size N. Fix an 
HVMS, {<y n } n€ j N i of size N + 1. The case when / = 1 is handled as in the 
previous paragraph. If 2 < / < N + 1, as {a; n } n6j is an HVMS of size N, 
the inductive hypothesis implies that 

Ri(z)a = Zy 1 ARi_i(z)a = Zy 1 A ^ ~^ a n- ( 2 -9) 

\n\=l-l 

As {<y n } ne j N j is of size N + 1 and I — 1 < N, 2.3 implies that a n G Dom(A) 
whenever | rz | = / — 1. Hence, 2.9 implies that a G Dom(i?/(z)). Also, using 
2.2 and 2.4 we see via 2.9 that 



a = z Y 1 A ^2 \ a n 

\n\=l-l 

= *r E ^o re+ei + ^ £ 1(1 -Y) 

|n|=Z— 1 |n|=i— 1 

1 ^ 



m|=i 



m^(Z,0),(0,Z) 

1 E ^(l-^+^ill-^M 

ii 2 2 
m^(J,0),(0,i) 



- r a (J)0 ) + z Y l Yl Z^ Ya ™+ 

1 \m\=l 

m^{l,0),(0,l) 



2 

|m|=Z " z 
m^(Z,0)(0,Z) 

n-a(i,o) + ^y 1 E (^i^ + ^(1 - Y))la m + -^-a (0 ,i) 



|m|=Z 
m^(Z,0)(0,Z) 



2. 



2 



s 




1 1 

— a m + -fa m 
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m^(Z,0)(0,Z) 



E 1 

[n|=J 



□ 



The property described by 2.7 in Proposition 2.5 arises as an issue in 
many of the applications of HVMS's that we have in mind. Accordingly, we 
introduce the following definition. 

Definition 2.10. Let H be a Hilbert space, a £ Ti, and assume that Y 
is a positive contraction on H. If A is a symmetric operator on H , we 
say that A has finite complex vector (Y, a) -moments to order N if for each 
z £ {z\z 2 ^ 0,zx/z 2 i (-oo,0]}, a £ Dom((A^ 1 ) / ) fori = 1,...,N. 
We say that A has finite real vector (Y, a) -moments to order N if for each 
b £ R+ 2 , a £ Dom((Aby 1 ) 1 ) for I = 1, . . . , N. 

The following converse to Proposition 2.5 provides a useful criterion to 
verify that a given symmetric operator and positive operator are associated 
with an HVMS. 

Proposition 2.11. Let H be a Hilbert space, let a £ H and assume that A 
and Y are operators acting on H, with A symmetric and Y a positive con- 
traction. The following conditions are equivalent. 

(i) There exists a sequence {« n } neJjv in H such that 
ol = tt (i,o) + a (o,i) an d ({ a n}, A, Y) is an HVMS. 

(ii) A has finite complex vector (Y, a) -moments to order N — 1 and 
for each I = 1, . . . , N there exist vectors a n , \n\= I such that 




\n\=l 



whenever z £ {z \ z 2 ^ 0, zi/z 2 £ (-co, 0]}. 
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(iii) A has finite real vector (Y, a) -moments to order N — 1 and 
for each I = 1, . . . , N there exist vectors a n , \n\= I such that 

whenever b G R +2 . 



Proof. That (i) implies (ii) follows from 2.5. Obviously, (ii) implies (iii). 

Assume that (iii) holds. To show that 2.2 holds when I = 1 and that 
ct = a (i,o) + Qf(o,i), equate coefficients in the following equation obtained from 
2.12 when 1 = 1. 

a = byby^ct 
= b Y Ri(b)a 

= {biY + 62(1 - K))(ia ( i, ) + r«(o,i))- 



% K ^ b 



2 



Now assume N > 2. Note that the moment condition implies that for 1 < 
I < N - 1, Ri(b)a G Dom(v4). Hence by 2.12, 



J^ a n e Dom(A), 



b 

\n\=l 



for all beR +2 . As 



span{^^ tyi an I ^ e R +2 } = span{a n | |n|= /}, 

|rc|=Z 

it follows that 2.3 holds. 

Now fix Z with 1 < I < N - 1. Noting that b Y Ri+i(b) = AR t (b), we 
compute using 2.12 that 

E- — Aa m = A - — a m 

\m\=l \m\=l 

= bY 2 ^ 

|n|=/+l 
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= ^T^"(0,m) + Yl ^ Y(Xm +^ + ( X ~ Y ) a m+e 2 ) 
2 \m\=l 

b 2 

+ ^Ti( 1 - y >(m,o)- 

Equating terms in this formula yields that 2.4 holds for 2 < I < N — 1 and 
that 2.2 for 2 < / < N. 

□ 

We now turn to a much more subtle characterization of HVMS's given 
in Theorem 2.21 below. Suppose that {a n } is an HVMS of size N with 
weight Y and shift A and let a be as in 2.6. Let R + = {t E R 1 1 > 0}. For 
1 < k < 2N — 1 define functions r k : IR +2 — > M by the formulas, 

ri(b) — < Ri(b)a, a > — < b^a, a > if k = 1 (2.13) 

r k (b) = < R t {b)a, ARi^{b)a > if 3 < k = 21 - 1 (2.14) 

r k (b) = < Ri{b)a, AR t (b)a > if 2 < A; = 21, (2.15) 

where the expressions Ri(b)a make sense by Proposition 2.5. Computing 
rfc(6) using 2.8 yields the qualitative information that for each k with 1 < 
k < 2N — 1, r k {b) is a homogenous polynomial in | = (^-,5-) of degree 
To formalize these properties of a, Y, and A we introduce the following 
definition. 

Definition 2.16. Let H be a Hilbert space, a G Ti, and assume that Y is 
a positive contraction on %. Assume that A is a symmetric operator on H 
with finite real vector (Y, a) -moments to order N — 1. For 1 < k < 2n — 1 
we define the k th scalar (Y,a)-moment of A by equations (2.13) to (2.15). 

Before continuing, we remark that ontologically the scalar (Y, a)-moments 
of A are functions on (M + ) 2 . However, if these functions happen to be given 
by homogenous polynomials (as e.g. occurs in the case of an HVMS), then 
there is an obvious way to extend the moment functions to all of C 2 . Concrete 
formulas for this case would be 

ri(z) = < R\(z)a,a > = < z^a^a > if k = 1 (2.17) 

r k {z) = <R l {z)a,AR l - 1 {z)*a> if 3 < k = 21 - 1 (2.18) 

r k (z) = < Ri(z)a, AR^zfa > if 2 < k = 21. (2.19) 
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Remark 2.20. If ({a n },A,Y) is a finite HVMS, then by Proposition 2.11 
the k th scalar (Y, a)-moments of A are given by 

r i( b ) = 67 ("(i,o), ") + ^(«(o,i),a) if fc = 1 

r k(b) = J2\m\=i-i,\ n \=i ^h^(a n ,Ya m+ei + (1 - Y)a m+e2 ) if 3 < k = 21 - 1 

r k(b) = J2\m\=l,\n\=l b^( a ^ Ya m+e 1 + (1 ~ Y)a m+e2 ) if 2 < fc = 21. 

In particular, they only depend on the Gram matrices a 1 = {Ya n , ct m ) and 
a 2 = ((1 - Y)a n ,a m ). 

Theorem 2.21. Let "H be a Hilbert space, a and N > 1. ^ss^me that 
Y is a positive contraction on % and A is a symmetric operator on H . There 
exists an indexed sequence {oc n } neI of vectors in H such that ({a n },A, Y) 
is an HVMS of size N and 

a = a ( i i0 ) + «(i,o) (2.22) 

if and only if A has finite real vector (Y, a) -moments to order N — 1 and 
for each k < 2N — 1, the k th scalar (Y, a) -moment of A is a homogeneous 
polynomial in | of order k. 

Proof. The necessity of the homogeneity condition follows by the discussion 
leading up to Definition 2.16. To prove the sufficiency we proceed by induc- 
tion on N. 

When N — 1, there is only one scalar moment given by 

ri(b) =< b Y l a, a > . 

If ri is homogenous of degree one, then there exist constants ai and a 2 such 
that 

< b^a, a >= ai— + a 2 — . (2.23) 
b\ b 2 

We analyze 2.23 by making the substitutions, 

bi = x and b 2 = — ^——x. (2.24) 

t 1 

Noting that in the new variables x and t, 

by^—it-Y)' 1 , (2.25) 
x 
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one computes that 2.23 becomes 



<(t-Y)- L a,a>= -^j + j. (2.26) 

Now, 2.26 implies that the scalar spectral measure of Y w.r.t. a is supported 
in the set {0, 1}, which in turn implies that Y(l — Y)a = 0. Letting «(i,o) — 
Ya and 01(0,1) = (1 ~~ Y)q:, we see immediately that 2.22 holds. As 

(1 - Y)a m = (1 - Y)Ya = 

and 

Ya m = Y(l - Y)a = 0, 

we see that 2.2 holds. Finally, as 2.3 and 2.4 are both vacuous when N — 1, 
the theorem is proved for the special case when N = 1 

Now suppose the sufficiency of the homogeneity conditions whenever A 
has finite vector (Y, a)-moments to order N—l and 2iV— 1 homogenous scalar 
(Y, a)-moments. Fix A, Y, and a with the properties that A has finite real 
vector (Y, a)-moments to order iV and 2N + 1 homogenous real scalar (Y, a)- 
moments, ruip). We need to show that there exists an indexed sequence 
{ a n} n€ j N+1 in T~L such that ({a n }, A, Y) is an HVMS of size N + 1 and such 
that a = «(i,o) + onm. By Proposition 2.11 this will be accomplished if we 
can construct an indexed set \a n }„ cT in H such that 

{a n } n£lN C Bom(A) (2.27) 

and 

R^a = J^ a n for / = 1, . . . , N + 1. (2.28) 

\n\=l 

By the induction hypothesis, there exists an indexed set of vectors in %, 
{a n } n( zj N such that 2.22 holds and such that ({o„,} ne/jv , A, Y) is an HVMS of 
size N. By the homogeneity of r 2 N+i(b), there exist scalars p n , \n\= 2N + 1, 
such that 

r 2N+1 (b)= J2 % ( 2 -29) 

|n|=2AT+l 

On the other hand, by the definition of the odd scalar moments, 2.14, 

r 2N+1 {b) =< by X AR N (b)a, AR N {b)a > . (2.30) 



13 



Finally, Proposition 2.5 implies that 

R N (b)a = ^«n- (2.31) 

|n[=JV 

The remainder of the proof consists of employing the substitutions, 2.24, to 
make various deductions from 2.29, 2.30, and 2.31 pertinent to establishing 
2.27 and 2.28. To facilitate our calculations we shall employ the notation, 

t n = t ni (t-l) ri2 . 
Making the substitutions, 2.24, in 2.31, we obtain that 

R N {b)a = [tx)~ N ( 2 - 32 ) 

\n\=N 

As A is assumed to have finite real vector (Y, a)-moments to order N, the left 
side of 2.32 is in the domain of A for all b G (1R + ) 2 . Hence, the right side of 
2.32 is in the domain of A for all t G 1R + . Noting that the set {t n \ \n\= N} is 
a basis for the polynomials of degree less than or equal to N, it follows that 
a n G Dom(A) whenever \n\= N. On the other hand, as ({a n } nelN , A, Y) is 
an HVMS of size N, it follows from 2.3 that a n G Dom(A) whenever \n\< N. 
Thus, we have shown that 2.27 holds. 

In order to verify 2.28 we must first explain how a n is defined when 
\n\= N + 1. Substitute 2.31 into 2.30 and then equate the right hand sides 
of 2.29 and 2.30 to obtain, 

(te) -<*m> £ pnt n 

|n|=2JV+l 

t-1 



<{t-YY\txY N t m Aa m ,{txy N r ^™> 

|m|=JV |m|=Af 



which simplifies to 

p(t) 



-< (t-Yy 1 t m Aa m , tmAa m >, (2-33) 

^ ' \m\=N \m\=N 

where p is the polynomial of degree less than or equal to 2N + 1 defined by 

p(t)= E p- tr ( 2 - 34 ) 

\n\=2N+l 
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For m a multi-index, we define Q m (t), an operator- valued polynomial of 
degree \m\ —1, by the formula, 

t m -Y m 

Qm\t) — 7 _ y ' 

Computing with the right side of 2.33 yields that 
<(t- Y)- 1 t m Aa m , tmAa m > 

\m\=N \m\=N 

= <(t- Y)- 1 K* - y)Qm(t) + Y m ]Aa m , tmAa m > 

\m\=N \m\=N 

= < J2 Qm(t)Aa m , > 
\m\=N \m\=N 

+ <(t- Y)' 1 Y m Aa m , t m Aa m > 

\m\=N \m\=N 

= < Y Qm(t)Aa m , t m Aa m > 

\m.\=N |m|=7V 

+ < Y ymA ^ (* - y y x Y tmAa m > 

\m\=N \m\=N 

= < Qm(t)Aa m , tmAa m > + < YmAQL ^ Y Qm(t) A a m > 

\m\=N |m|=7V |m|=7V \m\=N 

+ < Y Y m Aa m , (t - Y)- 1 Y Y m Aa m > . 

\m\=N \m\=N 

As the first two terms of this last expression are polynomials of degree less 
than or equal to 2N — 1 and N — l respectively, recalling that p has degree less 
than or equal to 2iV + 1, we see that the third term in the above expression 
must have the form, 

< Y Y m Aa m , (t - rr 1 Y YmAa m >=J + + ?(*)> ( 2 - 35 ) 

|m|=iV \m\=N 

where C\ and c<i are scalars and q is a polynomial of degree less than or equal 
to 2N — 1. 2.35 implies that if we set 

p= Y, y mA ^m 

\m\=N 
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and E is the spectral measure for Y, then dEp^ is supported in {0, 1}, which 
in turn implies that 

Y(Y-l)f3 = 0. (2.36) 
Now observe in light of 2.36, that 

t(t - 1)0 - IT 1 /? = (t + Y - l)/3 

Hence, 

i(*_l)(t_y)-i ^ t "M« m 

|m|=7V 

= ^-^(t-y)- 1 ]T [(t-F)Q m (t)+F m ]^«- 

= t(t - 1) ^ Q m (t)Aa m + t(t-l)(t-Yy 1 Y m Aa m 

\m\=N \m\=N 

= t(t-l) J2 Q m (t)Aa m + t(t-l)(t-YY 1 p 

\m\=N 

= t(t-l) £ Q m (t)Aa m + (t + Y -l)/3. 

\m\=N 

As Q m (t) has degree iV — 1, this implies that 

t(t- l)(t-y) _1 ^ t m Aa m 

|m|=JV 

is a vector valued polynomial of degree N + 1. But the set {t n | |n|= iV + 1} 
forms a basis for the polynomials of degree less than or equal to N+l. Hence, 
there exist vectors a n 6 H, \n\= N + 1, such that 

t{t- l)(t-y) _1 t m Aa m = ra "- ( 2 - 37 ) 

\m\=N \n\=N+l 

Unraveling the substitutions 2.24, and 2.37 becomes 

b^AR N {b)a = Y J^ a ^ 

\n\=N+l 

or, 

R N+1 (b)a = ( 2 - 38 ) 

\n\=N+l 
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In addition, recalling that ({a n } ne/jv , A, Y) is an HVMS of size N, we see 
from Proposition 2.5 that 

R t (b)a = J^ a n for I = 1, . . . , N. (2.39) 

\n\=l 

Taken together, 2.38 and 2.39 imply 2.28. □ 



3 From HVMSs to Loewner Functions 

We let IT 2 denote the set {z 6 C 2 Im(*i) > 0, lm(z 2 ) > 0}. We let V denote 
the Pick class on n 2 , i.e. the set of holomorphic functions on II 2 that have 
nonnegative imaginary part. If D C M 2 , we define the Loewner class, C(D), 
by 

C(D) = {h G V | h is analytic and real valued on D}. 

C(D), which captures a semi- local version of the notion of inner, arises in 
a variety of problems involving interpolation, the real edge of the wedge 
theorem, and the analysis of operator monotone functions — see e.g. [3]. In 
this section we wish to consider a fully local version of C(D). To that end 
we shall require a number of definitions. Let = Ik U {(0, 0)}. 

Definition 3.1. For x G M 2 and S C IT 2 let us agree to say that S approaches 
x non-tang entially , S ^> x, if x G S~ and there exists a constant c such that 

\\z — x\\ < c min{Im(2; 1 )Im(z 2 )} 

for all z G S 

Definition 3.2. Let h G V and x G 1R 2 . We say that x is a C k -point of h 
if h is "non-tang entially C k at x" i.e. there exists an indexed set of scalars, 
5 = {d~ n } n( - Jk , such that if S C II 2 and S —> x,then 

h(z) - V p T 5 n z n 
lim _J t=U*± h = o. (3.3) 

zes IFII 

Evidently, if h G C(D), x G D, and x is a C fc -point of h, then 5, as 
uniquely determined by 3.3, has the property that 5 n is real whenever n G 
Jfc. This suggests the following definition as a reasonable localization of the 
Loewner class. 
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Definition 3.4. Let k > and x G M 2 . I//iG V, we say that h is Loewner 
to order k at x if x is a C k -point of h and if 5, as uniquely determined by 
3.3, has the property that 5 n is real for all n e Jk- 

We introduce in Definition 3.9 below a class of functions, C N , obtained 
by adding three extra minor provisos to the notion in Definition 3.4. First we 
shall assume that k = 2N—1 is odd. Secondly, we wish to consider regularity 
as z approaches infinity non-tangentially rather than as z approaches a finite 
point x G IR 2 . Finally, we shall normalize h to have the value zero at infinity. 

To formalize regularity at oo, we introduce the following two definitions. 

Definition 3.5. If {z n } is a sequence in IT 2 , we say z n — > oo if z n = (X n , fi n ) 
and both X n — > oo and fi n — > oo. For S C n 2 we say that S approaches oo 
non-tangentially , S oo, if there is a sequence {z n } in S such that z n — > oo 
and a constant c such that 

\\z\\ < c min{Im(2; 1 )Im(z2)} (3-6) 

nt 

for all z G S . If S —> oo, we let adj(5) denote the smallest constant such 
that 3.6 holds for all z G S. 

Definition 3.7. Let Q be a metric space, u G Vt, and F : fl 2 — Q a map. 
We say 

F(z) lo as z oo 
if for each S C fl 2 such that S oo, 

lim F(z) = oj. 

zes 

We now can extend the notion of C fc -point to oo. 

Definition 3.8. If h £ V we say oo is a C k -point of h if there exists an 
indexed set of scalars, p = {p n } n£ j k , referred to as residues, such that 

\\z\\ k (h(z) — ^) —7-0 as z ^4 oo. 

neJ fe 

Finally, notice that the residue, P(o,o) ; when it exists, is the limit of h(z) 
as z — > oo non-tangentially, and hence we denote it by h(oo). Our third 
proviso is to normalize h by requiring that h(oo) = 0. 
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Definition 3.9. For N a positive integer, let £ denote the set of all h G V 
such that oo is a C 2N ~ 1 -point for h with real residues and h(oo) = 0. 

Let us note that Theorem 1.9 implies that any function in £} must have 
a representation as in (1.10). In previous work [3], we required Y to be a 
projection. This was inspired by representations on the bidisk, as in [6] and 
[2] . Here, we do not require Y to be a projection necessarily. But, in order for 
({a n } ng7iv , Y, A) to be an HVMS of size N, the operator Y{\ — Y) annihilates 
OL(ifi) and a(o,/) for 1 < I < N. So these vectors "think" Y is a projection. 

We now can formulate the main result of this section. 

Theorem 3.10. Suppose that H is a Hilbert space, A is a densely defined 
self-adjoint operator on %, a G %, and h G V is defined by the type I 
Nevanlinna representation, 

h(z) =< {A- z Y )~ 1 a,a>, zeU 2 . (3.11) 

If ({«„}, A, Y) is a of size N and a = acn^ + ccm.i), then h G C . Further- 
more, if Vj are the scalar (Y, a) -moments of A (as given by the formulas 2.18 
and 2.19) and p n are the residues of h (as given in Definition 3. 8), then 

(3 - 12) 

\n\=l 

fori = 1,...,2N - 1. 

The remainder of the section will be devoted to the proof of Theorem 
3.10. Accordingly, fix an HVMS of size N, ({a n },A, Y), with the property 
that A is densely defined and self-adjoint, set a = a(i,o) + a (o,i) ; an d assume 
that h is given by 3.11. The point z will always lie in n 2 , so (A — zy) is 
invert ible. 

Observe that as, 

A(A - z Y )~ l = (A-z Y + z Y )(A - zy)' 1 = 1 + z Y (A - z Y y\ (3.13) 

the operator A{A — z Y )~ l is bounded. Likewise, the operator (A — z Y ) _1 A 
is bounded. Also, we have the following simple identities involving these 
operators: 

z Y (A - z Y )~ l A = A(A - z Y )~ l z Y (3.14) 
(A - zyY 1 = -zy 1 + Zy X A{A - zyY 1 (3.15) 
(A - zyY 1 = -Zy 1 + (A - z Y y 1 Azy 1 . (3.16) 
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Claim 3.17. 

2N-1 

< (A — zy)~ X oc, a >= — rk{z)+ < A(A — z Y )~ l z Y R^(z)a, Rn(z)*oi > . 

k=l 

Note that as ({a n }, A, Y) is an HVMS, Condition (ii) of Proposition 2.11 
guarantees that a G Dom(R N (z)) and in addition, that the residues, r k (z), 
k = 1, . . . , 2iV — 1, are well defined by equations 2.17 to 2.19. Thus, the 
expression that appears on the right side of the claim is well defined. 

To prove Claim 3.17 we proceed by induction. Note that when N — 1 the 
claim follows immediately from 3.15. Suppose the claim holds for HVMSs of 
size N. If ({a n }, A, Y) is an HVMS of size N + l, then as ({a n }, A, Y) is 
also a HVMS of size N, the inductive hypothesis yields that 

2N-1 

< (A — 2y) -1 a!, a >= — rk(z)+ < A{A — zy)~ X zyRn(z)o!, RN{z)*a > . 

k=l 

(3.18) 

But, 

< A(A — zy)^ 1 zyRn(z)®, RN{z)*a > 
00 = < (A - z Y y 1 z Y R N (z)a, AR N (z)*a > 
(U) = < (-z' 1 + (A - z Y )~ 1 Az Y 1 )z Y RN{z)a, AR N (z)*a > 

= - < R N (z)a, AR N (z)*a > + < (A - z Y y 1 AR N (z)a, AR N (z)*a > 
(Hi) = — < Rn(z)ot, AR N (z)*a > 

+ < (-Zy 1 + z Y 1 A(A - z Y y 1 )AR N (z)a, AR N (z)*a > 
= - < R N (z)a,AR N (z)*a > - < Zy 1 AR N (z)a, AR N (z)*a > 

+ < z Y 1 A(A - z Y y l AR N ( z)a, ARn(z)*q > 
= - < R N (z)a,AR N (z)*a > - < Zy 1 AR N (z)a,AR N (z)*a > 
+ < A(A - z Y y 1 z Y Zy 1 AR N (z)a, Zy 1 *AR N (z)*a > 

(iv) = — < Rn(z)cx, ARN(z)*a > — < RN+i(z)a 1 ARn(z)* a > 

+ < A(A — z Y y 1 z Y R N+1 (z)a, R N+1 (z)*a > 

(v) = -t 2 n(z) - r 2 N+i(z) 

+ < A(A — z Y )~ z Y R N+1 (z)a, R N+1 (z)*a > . 
Here, the following facts were used. 
(i) as ({ctj, A, Y) is an HVMS of size N + l, R N (z)*a = R N (z)a e Dom(A) 
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(ii) 3.16 

(Hi) 3.15 

(if) RN+i(z)a = Zy 1 ARn(z) 

(v) 2.18 and 2.19 

Combining the result of this calculation with 3.18, we deduce that 

2JV+1 

< (A — z Y ) 1 a, a >= — 2_. Tk{z)+ < A(A — z Y )~ l z Y R N+1 (z)a, Rn+i(z)*ch >> 

fc=i 

which is 3.18 with N replaced with jV + 1. This concludes the proof of Claim 
3.17. 

Now observe that both the facts we need to prove to establish Theorem 
3.10, that h G C N and 3.12, will follow from Claim 3.17 if we can show that 

\\z\\ 2N ~ l < A(A- z Y y 1 z Y RN(z)a,R N (z)*a>^$ as z 4 oo. (3.19) 
On the other hand, we claim that 3.19 will follow if we can show 
Claim 3.20. 7 G U, then 

< A(A - z Y y l /3, 7 as z4oo. (3.21) 

To see how Claim 3.20 implies 3.19 we use the following simple property 
of sets that approach 00 non-tangentially. 

Lemma 3.22. If n is a multi-index, S C II 2 and S 4 oo ; then 

|l|<(adj(^))H|| z ||-N 

for all z G S . 

Proof. If z G then Definition 3.5 implies that 

IMI < a dj(5) min{Im(2; 1 )Im(^2)} 

< adj(S') Im(,2i) 

< adj(S') \z x I . 

Hence, 

|*i I"" 1 < adj(S)"' 1 \\z\\- ni . 

Likewise, 

\z2\~ n2 < adj(S)" 2 ||2|r n2 . 
The lemma follows by multiplying these last two inequalities together. □ 
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Now, using Proposition 2.5, if iV > 2, 

< A(A — z Y )~ zyRn{z)oi-, R,N{z)*ct > 
= < A(A - z Y y 1 AR N _i(z)a, R N (z)*a > 

= <a(a- z Y y x A Y ~^ a m, Y i a « > 



| m |=JV-l \n\=N 



Y -^ <A ^ A ~ z y) X Aa m ,a n >. 



z 

\m\=N-\ 
\n\=N 

Thus, using Lemma 3.22, we see that if S —)■ oo and z G S, then 

| < A(A — zy)' 1 zyRn{z)c(, RN{z)*a > \ 
^ Y I <A(A-z Y y 1 Aa m ,a n > \ 



\m\=N-l 
\n\=N 



< adj^) 2 ^ 1 !^!!-^- 1 ) Y \<MA-z Y y 1 Aa m ,a n >\. 



|m|=JV-l 
|n|=JV 



When N — 1, we get 

| < A(A - z Y y 1 z Y Ri(z)a, Ri(z)*a > 

= i Y ~^ < A ( A ~ ^r 1 "' ° n > i 



\n\ = l 



< ^{s)\\ z \y l Y\ <A ^ A 

- z Y ) a, a n >\. 



n =1 



So we see that Claim 3.20 does indeed imply 3.19. 

There remains to prove Claim 3.20. For this we shall require three lem- 
mas. These lemmas involve the notion of a proximity estimate, an idea which 
we make precise in the following definition. 

Definition 3.23. Let Q be a metric space and F : LT 2 — > Q a map. We say 
that F is proximal (or more precisely, proximal at oo) if for each S C n 2 

nt 

such that S — > oo, there exists a constant c such that 

d(F(z), F{w)) < c " p (3.24) 
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for all z,w G S. We refer to the inequality 3.24 as a proximity estimate. 

It turns out that frequently, as a consequence of various forms of the 
Schwarz Lemma, quantities that are formed from holomorphic functions sat- 
isfy proximity estimates. In such cases, the following lemma greatly simplifies 
the analysis of non-tangential regularity. 

Lemma 3.25. Let Q be a metric space, let to G Q and let F : II 2 — > Q be a 



proximal map. F(z) — > to as z ^> oo if and only if for each 5 G II 2 7 

lim F(s5) = to. (3.26) 

s— too 

Proof. Clearly, if F(z) — > to as z oo, then 3.26 holds. To prove the 
converse we argue by contradiction. Suppose 3.26 holds. If it is false that 
F(z) w as 2 4 oo, then there exist e > 0, S C II 2 , and a sequence {z{\ in 
S such that S oo, z\ — > oo, and 

d{F{zi),u)>e (3.27) 

for all positive I. By compactness, there exist 5 G C 2 and a subsequence 
zj., such that \\zi d H -1 ^ —> 8 as j — > oo. In fact, 5 G II 2 . To see this, let 
5 = (Si,8 2 ) and z\ } = (Aj.,^) and observe that 

Mb) = 1- > 1- min{Im M! mK)} > > 0, (3.28) 

3-*°° \\zijW i->°o ||^. || adj(6) 

Likewise, 1100.(82) > and we conclude that 5 G II 2 . 

Now, let Wj = zij and Sj = \\z^ ||. By construction we have that Wj — Sj5 = 
o(sj) so that the proximity estimate gives that 

Win ■ — s -i^ 

d(F(wj), F( Sj 8)) < c 11 3 3 11 ->• 0. (3.29) 

|| SjO || 

Also, 3.26 implies that 

d(F(sj-5),w) ->-0. (3.30) 
Given 3.29 and 3.30, the triangle inequality gives that 

d(F( Wj ),u))^0, 

contradicting 3.27. □ 
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Lemma 3.31. Let C(H) denote the algebra of bounded operators on % 
equipped with the operator norm. F : U 2 — > C(H), defined by 

F(z) = A(A — z Y y l ,zeU 2 , 

is proximal. 

Proof. Fix Sell 2 with S -^f oo. For z G II 2 we have that 

\\z Y \\ < max{|zi|, \z 2 \} < y/2\\z\\. (3.32) 

Also, as Im(v4 — Zy) = —hn(zy) < — min{Im(zi)Im(£2)}, we have that if in 
addition, z G S, then 

II {A - zy)- 1 1| < . n ( \. . u < (3.33) 
mm{Im(z 1 )Im(z 2 )} |p|| 

Now, using (3.13), we get 
- F(w) 

= A(A - z Y y l - A(i4 - wy) -1 

= (1 + Zy(A - Zy)' 1 ) - (1 + Wy(-4 - Wy)" 1 ) 
= Zy(A — Zy)' 1 — Wy(A — Wy)" 1 

= (Zy - Wy)(A - Zy)~ + Wy((A - Zy)^ - (A - Wy) _1 ) 

= (zy — Wy)(A — Zy)~ l + Wy(A — Wy)~ (zy — Wy){A — Zy)^ 1 . 

Hence using 3.32 and 3.33, 

\\F(z) - F(w)\\ 

< \\zy — Wy\\ \\(A — Zy)~ \\ + \\Wy\\ \\ (A — Wy)~ l \\ \\zy — Wy\\ \\ (A — Zy )~ 1 

^ /on ii ad K g ) , /on ii ad K^) /on I, adj(S0 

< V2 2 — u> — — h V 2 ky — — — V2 2 — u> — — — 

\\z\\ \\w\\ \\z\\ 



(V2 adj(5) + 2 adj(5) 



2\ 



\Z — W\ 



which is 3.24 with c = v / 2adj(5) + 2adj(5) 2 . □ 
Lemma 3.34. I//3,7 G H and 5 G II 2 , then 

lim < v4(A - t5 Y y 1 ^ 7 >= 



24 



Proof. We claim that for each vector u & W, 



Sy(eA — 8y) u — > —u weakly in V, (3.35) 

as e — > 0. To prove this claim first notice that as lm(5y) > min{Im(5i)Im((52)}, 
we have both that Sy is invertible and that (eA — 5y)~ 1 is uniformly bounded. 
In particular, as A is densely defined, M. = 5yDom(A) is dense in "H. If 
u = 5yv G Ai, then as v G Dom(A) and (eA — 5y)~ 1 is uniformly bounded, 



(eA - Sy^eAv 



as e — )• 0. Hence, 



(eA - = (eA - 5y)- l 5 Y v 

= {eA - 5 Y y 1 {{5Y - eA) + eA)v 
= -v + (eA - 5 Y y 1 eAv 
-> — v 
= —Sy^u. 

Applying the bounded operator 5y, yields that 3.35 holds whenever u G M.. 
As, M. is dense and 5y(eA — <5y) _1 is uniformly bounded, it follows that 3.35 
holds for all u G H. This proves the claim. 

Now notice that if in the claim, we substitute e = s -1 , we deduce that 
for all u 6H, 

sSy(A — s5y)~ 1 u — > —u weakly in 7i 
as s — > oo. Hence, for all «6H, 

1 + s5y(A — s<5y) _1, u weakly in "H 

as s — >■ oo. The lemma now follows by observing that from (3.13) 

1 + s5y(A - Sdy)- 1 = A(A - Sdyy 1 

□ 

Armed with the above lemmas it is a simple matter to prove Claim 3.20 
and thereby complete the proof of Theorem 3.10. If /3,7 G H, then by 
Lemma 3.31 F(z) =< A(A — zy)~ l f3,^f > is proximal. As Lemma 3.34 gives 
that lim^oo F(s5) = whenever 5 G n 2 , Lemma 3.25 yields that F(z) — > 
as z oo as was to be proved. 
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4 From Loewner Functions to HVMSs 



In this section we shall formulate and then prove a converse to Theorem 3.10, 
using Theorem 2.21. If h G C N , then it is easy to check that h is type I and 
accordingly has a Nevanlinna representation of the form, 

h(z) =< (A - z Y y l at, a >, z G II 2 , (4.1) 

where A and Y are operators acting on a Hilbert space H, A is a densely 
defined and self-adjoint, Y is a positive contraction, and a£H. 

Theorem 4.2. If h £ C N and A, Y , and a are such that 4-1 holds, then 
A has real vector (Y, a) -moments to order N — 1 and homogenous scalar 
(Y, a) -moments to order 2N — 1. Furthermore, 

whenever 1 < I < 2N — 1 and 6 G (ffi + ) 2 , where p n are the residues of h. 

Proof. We proceed by induction. Let N = 1 and assume that h G has 
a Nevanlinna representation as in 4.1. As N = 1, the assertion that A have 
real vector (Y, a)-moments to order iV — 1 is vacuous. To see that A has 
homogenous scalar (Y, a)-moments to order 2N — 1, first note that since oo 
is a C^-point for h with real residues, we have that there exist P(i,o), P{o,i) £ ^ 2 
such that 

< (A-z Y y 1 aia >=h(z) = ^ + ^ + o(\\z\\- 1 ), (4.4) 

Zi z 2 

non-tangentially at oo. Fixing b G (R + ) 2 and setting z = isb in 4.4 gives that 
is < (A - i sb Y y 1 a, a > + (4.5) 

as s — 7- oo in IR + . Noting that for 6 G (M + ) 2 , 6y is strictly positive definite 
and hence, invertible, we define a self-adjoint operator, Xj,, by the formula, 

X b = b Y ^Ab Y K (4.6) 

Noting that 

-l 1 1-1 

is < {A — isby) a,a > = is < (b Y (Xb — is)b Y ) a,a > 
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is < b Y 2 (Xf, — is) b Y 2 a,a> 

is < (Xb — is) b Y 2 a, b Y 2 a > 
X b +is ,-i 



«s < 



X 2 + s 2-y 



-s 2 + isX bl -i 



XI + s 2 



o y 2 a, by 2 a >, 



we see upon taking real parts in 4.5 that 

< X 2 + s2^ "'^ h + b 2 [ {) 



as s — > oo in M + . Now, the Lesbesgue Dominated Convergence Theorem 
guarantees that 



s 2 



*6 +* 2 

as s — >■ oo in R + . Hence, 



6 v 2 a — >■ 6v 2 a 



i 



- < b?a, a>=-< by* a, b^a >= ^ + (4i 

As 4.8 holds for all 6 G (M + ) 2 , we conclude that A has homogenous scalar 
(Y, a)-moments to order 1 as was to be shown. Also note that 4.8 implies 
that 4.3 holds. 

We now turn to the inductive step of the proof. Accordingly, assume that 

A has real vector (Y, a)-moments to order N — 1, (4-9) 
A has homogenous scalar (Y, a)-moments to order 2N — 1, and (4.10) 

1 < I < 2N - 1, 6 g (R + f => % = ~ ri ^ ( 4n ) 

\n\=l 

whenever h G C N and has a representation as in 4.1. Fix h with a represen- 
tation as in 4.1 and assume that h G C N+1 . We need to show that 4.9, 4.10, 
and 4.11 hold with N replaced with N + 1. However, as h G C N+1 C C N , the 
inductive hypothesis implies that 4.9, 4.10, and 4.11 hold for N. Therefore, 
the induction will be complete if we can show the following three conditions: 

a G Dom((A6 y 1 ) iV ), (4.12) 



27 



T2N + iib) = - & e(M + ) 2 , (4.13) 



|n|=2JV+l 

and 

r2N{b) = " E F' bE (R+)2 ' (414) 

|ra|=2iV 

First note that as h G C N+1 and 4.1 holds, there exist scalar residues, p n , 
n G I2N+1, such that 

<(A-^)- 1 a,a>= £ ^ + o(||,|r (2Af+1) ) (4.15) 

"S/2iV+l 

as z — > 00 non-tangentially in II 2 . Fixing b G (M + ) 2 and setting z = isb in 
4.15 we deduce that 



2N+1 

<(A- i sb Y y 1 a, a>=J2 (is)' 1 £ £ + ^- (2iV+1) ), 



Z=l |n|=i 

as s — )■ 00 in IR + , which, upon taking the imaginary parts, yields that 

N+i /_ 1 a 

im(<(A-^ y rv«>) = E^r E ^+^~ (2Ar+1) ) ( 4 - 16 ) 

fc=l |n|=2fc-l 

as s 00 in R + . Finally, upon multiplying 4.16 by the factor s 2N+l , we 
deduce the limit, 

lim G b (s) = (-lf +1 E ^ (4.17) 

|n|=2AT+l 

where for s G M + and b G (M + ) 2 , G&(s) is defined by 

N 

G b (s) = s 2N+1 lm(<(A~ i sb Y y 1 a,a>)~Y^(- 1 ) ks2iN ~ k+1) E ^ 

fc=l |n|=2fc-l 

(4.18) 

We now compute G&(s) using the substitution 4.6. We set 

lh = b'^a. (4.19) 
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Note that 4.9 implies that 7 b G Dom(X^) for I = 1, . . . , N - 1. Using 4.11 



and 2.14 we see for k = 1, . . . , N, that 

|n|=2fe-l 

= - < by 1 {Ab Y l ) k ' l a, {Aby^f^a > 

= -< (x b ) k -\, (x b ) k -\ > . 

Also, just as in the calculation leading up to 4.7 we compute that 

t a ■ i. \-i X b + is -i 

< {A — tsby) a,a>=<— 2 -b Y ^a,b Y a>, 

X h + s 



Im(< (A — isby) «,«>)=<— 2 -b Y 2 a i b Y 2 a > 



so that 

-i . , _ s 

^2 

=< ^T^T6,76>. 

Hence we have that 

Q 2N+2 N 

We claim that the above sum telescopes. Indeed, using the fact that 
< (X b ) k ~\, {X h ) k ~ l lh > 

we compute that 

2iV+2 ^ 

G*00 = < ^7,, 7, > +^(-l)V^" fc+1 ) < (X 6 )*- l 7h, (X*)*- 1 ^ > 

6 fc=i 

s 27V+2 

= <xJT7^ b> 
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s 2N (< 



X 2 b 

xJT7^ b> + < 



+ s 2N ~ 2 (< 



X 2 + S 2 X ^ X ^ > 



< 



XI + s> 



X b j b ,X b j b >) 



Y 2 

i / -i\N 2/^ ^b / V \N-1 /v\N-l 

+ (-1) s (< - (Xfe) 7b , (X 6 ) 7 b > 

^b T s 



+ < y2 i g2 yXb) lb, (X b ) 7fe >) 

b 

Y 2 

' i\N 2 „ ^b /v\N-l /v\N-l 
~^ S < Y 2 I .2 76, M 7b >• 



This last calculation makes sense since 



XI + s 2 



is a bounded operator and 7b G Dom(AV) for I — 1, . . . , iV — 1. 

Now recall 4.17. From the formula for G b {s) just derived, we see that 

Js<3^W" 1 *W" I '» > =- S f- (420) 

6 |n|=2JV+l 

As X b is self-adjoint and 7b G Dompf 6 we can apply the spectral theo- 
rem to X b and thereby obtain the scalar spectral measure of 7b, \i. Analyzing 
the very existence of the limit on the left side of 4.20 in the space L 2 (fi) yields 
via the Lesbesgue Dominated Convergence Theorem that 

7b G Dom(X b N ). (4.21) 

Unraveling 4.21 via 4.6 and 4.19 gives that, 

a G Dom((A6 y 1 ) Af ), 

which is 4.12. Note also from 4.20 we have that 

< (X b ) (Xb) 7b, (X b ) 7 b >= - 2^ ^T> 

|n|=2JV+l 
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which unravels to 

w&) = - e (422) 

|n|=2AT+l 

which is 4.13. 

There remains to check 4.14. This is done by following the same line of 
reasoning that led from 4.15 to 4.20. One starts with 4.15 but with 2N + 1 
replaced with 2N: 

<(A-z Y y 1 a,a>= J2 ^ + o(\\z\r (2N) ) (4.23) 

ne/ 2 JV 

Proceeding as before, for a fixed b G (IR + ) 2 and s G 1R + we set z = isb in 
4.23. However unlike before, where we took imaginary parts to obtain 4.16, 
we now take real parts. This results in 

N /_-,\k 

Re(< (A - zsby)- 1 ^ a>) = J2 [ -^rJ2^ + °( s '^ ( 424 ) 

k=l \n\=2k 

as s — > oo in R + . Finally, upon multiplying 4.23 by the factor s 2N (rather 
than s 2N+1 as before), we deduce the limit, 

lim F b (s) = {-If £ ( 4 - 25 ) 

\n\=2N 

where for s G M + and b G (M + ) 2 , F b (s) is defined by 

JV-l 

F b (s) = s 2N Re(< {A-isbyy^co) - ^ (-1) ^ ^. (4.26) 

k=l \n\=2k 

Carrying out the telescoping argument, one computes that 

F 6 (s) = (-1) S < 2 {X b ) 76, (X 6 ) 7fc >, 

which implies via 4.25 the existence of the limit 

6 |n|=27V 
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As 4.21 holds, 4.27 implies that 

< (x b f(x b f-\, (x b f-\ b 

\n\=2N 

As this last equation unravels via 4.6 and 4.19 to 

r 2 N(b) = - ^' 

|n|=27V 

the proof that 4.14 holds is complete. □ 

5 Finite Hankel Pairs 

In this section we give an alternate matrix theoretic treatment of HVMS's 
based on the fact that it is possible to cleanly characterize the Gram matrix 
formed from the moment vectors of an HVMS. 

For X a set, we let £ 2 (X) denote the Hilbert space of square summable 
complex valued functions on X. If / G £ 2 (X), we let supp(/), the support of 
f, denote the subset of X defined by 

supp(/) = {x G X | f(x) + 0}. 

By a matrix on X we mean a square array of scalars, doubly indexed by 
the elements of X. If a = [a XjV ] is a matrix on X, then a induces a densely 
defined linear operator, also denoted by a, on the finitely supported functions 
in £ 2 (X) by the formula 

( a f)( x ) = a *,yf(y)- 

J/Gsupp(/) 

If a = [a X) y] is a matrix on X, then we say that a is symmetric if 

a>x,y = a-y,x for all x, y G X, 

and we say that a is positive semi-definite if for each (finite) choice of ele- 
ments, Xi,x 2 , ■ ■ ■ ,xi G X, and each choice of scalars, ci, c 2 , . . . , q G C, 

i 



32 



In this section we shall be exclusively interested in the case where X = 1^, 
for iV a positive integer. Note that naturally, if M < N, then £ 2 (Im) Q 
£ 2 (In), and in addition, that there is a pair of shift operators, S\, S2 '■ 
l 2 {I N -i) -)• £ 2 (I N ) defined by 

(Sxf)(n) 

(S*f)(n) 

If ({a n }, A, Y) is an HVMS of size N, then we may define a pair of matrices 
a = (a 1 , a 2 ) on I N by 

a m,n =< Ya n , a m > and a 2 m n =< (1 — K)a n , a m > for m,n E In- (5.1) 

Definition 5.2. M^e say that a = (a 1 , a 2 ) is a finite Hankel pair of size N if 
a 1 and a 2 are matrices on 1^ and there exists an HVMS of size N such that 
(5.1) holds. 

In Theorem 5.7, we give a characterization of when a pair of matrices is a 
finite Hankel pair. To see how this is a two variable version of Theorem 1.7, 
let us restate that theorem more abstractly. Let S : £ 2 ({0, 1, . . . , N — 2}) — > 
£ 2 ({0, 1, . . . , N - 1}) be the shift defined by Sf(j) = f(j - 1), j > 0, and 
Sf(0) = 0. 

Theorem 5.3. Let H be an N-by-N matrix. There is a self-adjoint operator 
A and a vector a with a G Dom(A k ) for 1 < k < N — 1 such that 

Hij = < A j a, A { a> < i, j < N - 1 

if and only if the following three conditions obtain. 

H is positive semi- definite. (5-4) 

H i+1J = H id+1 0<i,j<N-2. (5.5) 
Bupp(/) G {0, . . . , iV - 2} and Hf = => HSf = 0. (5.6) 

Here is our two variable version of Hamburger's Theorem 1.7. 
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J /(n- ei ) 

I 

f /(n-e 2 ) 

I 



n - e x G Jiv-i 
else 

n - e 2 G Jat-i 
else 



Theorem 5.7. Let a be a pair of matrices on 1^. Then a is a finite Hankel 
pair of size N if and only if the following four conditions obtain. 

a 1 and a 2 are positive semi- definite. (5.8) 

<4+ei,n + a m+e 2 ,n = «m,n+ ei + a m,n+e 2 whenever m, 71 G IjV-1- (5.9) 

a ( ,o,(o,o = a lmm = ^ or z = x > • • • ' N - ( 5 - 10 ) 

supp(/) G Jjv-i and (a 1 + a 2 )f = =4> (a 1 ^ + a 2 S 2 )/ = 0. (5.11) 

Proof. (Necessity) Assume that F, A) is an HVMS and (5.1) holds. 

Then (5.8) holds because F and 1 — Y are positive operators. (5.9) holds 
because the left-hand side is 

by (2.4). But the right-hand side of (5.9) by a similar calculation is (Aa n , a m ), 
which is equal to (a n , Aa m ) because A is self-adjoint and a m , a n are in its 
domain for m, n G ijv-i- Condition (5.10) follows from (2.2). 
Finally, if supp(/) G Ijv-i an d (a 1 + a 2 )f = 0, this says that 

( f( n ) a n,o:m) = 

for all m G /at. But 

(a 1 ^ + a 2 S 2 )f(m) = ^ f{n){Ya n+ei + (l-Y)a n+e2 ,a m ) 

nG/jv-i 

= ( A ( Yl f( n ) a n),a m ) 

= o, 

so (5.11) holds. 

(Sufficiency). Assume (5.8) — (5.11) hold. Choose vectors a n in a Hilbert 
space H so that their Grammian equals the matrix a 1 + a 2 : 



i \ l 



2 

m,n' 



Since a 1 < (a n , a m ), there is a positive operator F satisfying (5.1). Equation 
(2.2) follows from (5.10). 
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If N = 1, we can define A arbitrarily, e.g. by A = 0. 
If N > 2, we define A on the span of {a n }nei N -i by 

Act„ = y« n+ei + (1 - Y)a n+e2 . 

To check that this is a well-defined linear operator, we need to know that if 

/ J c n a n = 0, 

then 

} u Cn(Ya n+ei + (1 - Y)a n+e2 ) = 0. 
This follows from (5.11). It follows from (5.9) that A is symmetric. □ 



6 Infinite sequences 

As in one variable, passage from the finite to the infinite case is straight- 
forward and leads to some simplifications. Let I denote the set of pairs of 
non-negative integers, excluding (0,0). 

Definition 6.1. An infinite Hankel vector moment sequence is a 3-tuple, 
({«n}iiei, Y, A) where: {a n } n£ i is a sequence of vectors in some Hilbert space 
%; Y is a positive contraction acting on 7-L, satisfying for each I > 1 

Ya(o,i) = = (1 - Y)a { i fl) = 0; 

A is a densely defined self-adjoint operator on H with the property that 

{a n | n e 1} C Dom(A); 

for each n G X, 

Aa n = Ya n+ei + (1 - Y)a n+e2 . 

Theorem 1.15 becomes a description of functions in C°°. 

Theorem 6.2. A Pick function h of two variables has an asymptotic expan- 
sion 
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as z — > oo, for some real numbers p n , if and only if it has a representation 
as in (1.10) and for every such representation there is an infinite HVMS 
({«n}neX) A) with a = a(i,o) + a (o,i)- Moreover, {p n } are given by 

\n\=l 

whenever I > 1 and b G (M + ) 2 . 

Sufficiency of the condition follows from Theorem 3.10; necessity follows 
from the constructive proof of Theorem 4.2. 
We define an infinite Hankel pair by 

Definition 6.3. We say that a = (a 1 , a 2 ) is a infinite Hankel pair if a 1 and 
a 2 are matrices on X and there exists an infinite HVMS such that 

a ln,n = (Y®n, «m) and a^ n = ((1 - Y)a n , a m ) for m,n el. (6.4) 

If (5.9) holds for all N, then (5.11) holds automatically. So the infinite 
Hamburger theorem becomes 

Theorem 6.5. Let a be a pair of matrices onX. Then a is an infinite Hankel 
pair if and only if the following three conditions obtain. 

a 1 and a 2 are positive semi- definite. 

a m+e 1 ,n + a m+e 2 ,n = ° J m,n+e 1 + a m,n+e 2 whenever 711, U E X. 

a (o,o,(o,o = a («,o),(z,o) = f or I — 1- 
Here is a two variable version of Kronecker's theorem. 

Theorem 6.6. Let h 6 C°° . Then there is an infinite HVMS ({a n } ne j, Y, A) 
with a = ct(i.o) + 0(0,1), satisfying rank(a 1 + a 2 ) < oo and 

h(z) = ((A — zy)~ 1 ol, a) (6.7) 

if and only if h is a rational function. 
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Proof. If h is rational of degree (d\, 62), then by Theorem 8.12 h has a repre- 
sentation (6.7) on a Hilbert space H of dimension at most d = d\ + di- Since 
h G C°°, by Theorem 6.2 there is an infinite HVMS ({a n }nex, Y, A) on TL. 
So 

(a 1 + a 2 ) mi „ = (a n ,a m ) n 

has rank at most 0?. 

Conversely, suppose there is an infinite HVMS ({a n }nei, Y, A) with a = 
a (i,o) + ^(0,1) j satisfying rank(a 1 + a 2 ) = d < 00 and (6.7). Then one can 
choose vectors (3 n in a space H of dimension <i such that 

(Pn,Pm)n = («n,a m ) (6.8) 

and so that the vectors {f3 n } span Ti. Define a positive contraction X on % 
by 

(Xp n ,p m ) = (Ya n ,a m ). (6.9) 

Define i? by 

-B/3 n = X(3 n+ei + (1 — X)(3 n+e2 . 

We claim that i? extends by linearity to a well-defined linear operator on %. 
Indeed, suppose J2 c n(3 n = 0. Then by (6.8), 

So 5^c n a n = 0, and therefore by (6.9) 

(^2c n [X/3 n+ei + (1 - X)p n+e2 ],p m ) = (y] c n [Ya n+ei + (1 - y)q n+e2 ], a m ) 

{A ^ ^ C n O n , Oirri) 

= 0. 

As {/3 m } span 7-L, this means ^c n [X/3 n+ei + (1 — X)(3 n+e2 ] = 0; so B is 
well-defined, and hence ({/3 n } ng x, X, £?) is an infinite HVMS on H. Let 

/3 = P(i,o) + /3(o,i)- 

By Remark 2.20 the scalar (X,/3) moments of B agree with the scalar 
(Y, a) moments of A to all orders. Therefore by Theorem 6.2, the rational 
function g of degree at most d in each variable given by 

g(z) = ((B-z x )- 1 (3,f3) n (6.10) 

has the same asymptotic expansion at oo as i By Lemma 6.11, we are 
done. □ 
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Lemma 6.11. Let g,h be in C°° and have the same asymptotic expansion at 
oo. Assume in addition that g is rational. Then g and h are equal. 

Proof. For each fixed w in K, the functions g(z, z + w) and h(z, z + w) are 
in the one variable Pick class and have the same asymptotic expansions at 
oo. By Theorem 1.4, they must be Cauchy transforms of measures with the 
same moments. Moreover, g(z, z + w) is rational. Therefore by [17, Thm. 
1.2], the one- variable moment problem is in this case determinate, so the 
two measures must be equal. Therefore g(z, z + w) = h(z, z + w) for all 
z G II, w G M, and so the two functions are identically equal. □ 

Corollary 6.12. Let h G C N have an asymptotic expansion 



as z — > oo. Then there is a rational function g in C°° that has the same 
asymptotic expansion to order 2N — 1 . 

Proof. Let {{c^n} n( zj N , Y, A) be a finite HVMS corresponding to h as in The- 
orem 4.2. Choose vectors {(3 n } n( zi N in a finite dimensional space % so that 
(6.8) holds, and define X and B as in the proof of Theorem 6.6. Then g 
given by (6.10) has the same asymptotic expansion. □ 

7 An example 

Let {vjj}'jL 1 be a summable sequence of non- negative numbers, and let {Xj}^ =1 
be a sequence of real numbers. Let tj be numbers in the interval [0, 1]. Define 




n|<2AT-l 



nt 
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"(0,1) 



e 



If a = a(i,o) + «(o,i) and h(z) = ((A — zy) a, a), then h(z) is given by 



in 

- J A| - A,(2^ - - z 2 ) - Z!^ 2 ' 



At^/l-t 2 X j + z 1 + z 2 



If ^WjX 2 ; < oo, then one can extend the HVMS by 



/ t* \ 



"(2,0) 



0(1,1) 



©v^^ (2^-1) 



1 - t 2 



02^ A, 



«(o,2) = Xj (1 - 2t)) 
Calculating, one gets that 
n{z) -- 



1-t 2 



r 2 {z) 



r 3 (z) 



1 1 

zi z 2 



2*3 - 1 



ZlZ 2 



3 1-2% 

z l 2 



(7.1) 



(2g - l) 2 4(^-^ + 4^-1)^1-^ 

z{z 2 



4(t 2 -t*)+4 tj (l-2t 2 )^l-t 2 (2 ^_i) S 

+ 5 + " ' 



^1^2 



7 3 
-2 



These are (up to a minus sign) the first 3 terms in the asymptotic expansion 
of (7.1) at infinity. If one assumes that Yl w j^j < °°> then one gets two more 
terms, and so on. 
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In the special case that every tj = 1/ a/2, the formulas simplify. Then 

k{z) = \)-z x z 2 



r 2 (*) = (j^A' 
^3(2) = (X) w J' A i 



2 



zi 22 
gi + z 2 
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8 Models 



A model for /i is a reproducing kernel space Ai on II 2 , and a positive con- 
traction F on Ai so that, if the reproducing kernel K for Ai is written as 

K(z,w) = (v z ,v w ) M (8.1) 

with v z analytic in z, then 



h(z)-h(w) = (zi-w 1 )(Yv Z: v w ) + (z 2 -w 2 )((I -Y)v z ,v w ). (8.2) 
Using our earlier notation zy = z{Y + z 2 (I — Y) , (8.2) becomes 



h(z)-h(w) = ((z Y - Wy)v z ,v w ). (8.3) 

The existence of models for functions in the Pick class was proved in [1]. 
Indeed, it was shown there that for every h in the Pick class, there are 
analytic functions v 1 (z) and v 2 (z) taking values in Hilbert spaces Ai l and 
M 2 so that 



h(z)-h(w) = (z 1 -w 1 )(v 1 (z),v 1 (w)} M i + (z 2 -w 2 )(v 2 (z),v 2 (w)) M 2. 
Let 

K(z,w) = (v 1 (z),v 1 (w)) M i + (v 2 (z),v 2 (w)) M 2. 

This is a kernel, so can be written as in (8.1) for some other Hilbert space 
Ai, and there is a positive contraction Y on Ai so that 

(v 1 (z),v 1 (w)) M i = (Yv z ,v w ). 

This yields (8.3). 

Write i for the point in C 2 . The equivalence of (ii) - (iv) in the 
following theorem was first proved in [4]. 
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Theorem 8.4. Let h : II 2 — > U be in the Pick class, and not identically zero. 
The following are equivalent. 

(i) For some/every model with reproducing kernel as in (8.1), there is a 
vector a in M such that 

h{z) = (v z ,a). (8.5) 

(ii) There exists a self-adjoint operator A on a Hilbert space H and a 
vector a in% such that 

h(z) = ((A - z Y )~ 1 ct, a). (8.6) 
(Hi) There exists c > such that 

lim sh(si) = ic. (8.7) 

s— >oo 

(iv) We have 

liminf \sh(si)\ < i. (8.8) 

Proof, (i) (ii): Define B by 

B : v z !->■ zyv z + a. 

Equations (8.3) and (8.5) imply that 

(Bv z ,v w ) = (v z ,Bv w ). (8.9) 

Extend B to finite linear combinations of vectors v z . by linearity, and (8.9) 
says that B is well-defined and symmetric. Indeed, if some linear combination 
c j v Zj = 0) then for every w we have 

(^2cj((zj) Y v Z:j +a),v w ) = (^2cjV Zj ,w Y v w + a) = 0, 

so B(cjV z .) = 0. 

If the defect indices of the closure of B match, then B can be extended 
to a self-adjoint operator on A^. If not, B can be extended to a self-adjoint 
operator on a superspace of Ai. In either event, we can assume that there is 
a self-adjoint A on "H 3 M. such that 

A:v x \-t z Y v z + a. 

Therefore v z = (A — zy)~ 1 ct, and (8.6) follows from (8.5). 
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(ii) =>- (Hi): By the spectral theorem, 

sh(si) = J —dfx(t) 

where \x is the finite measure that is the scalar spectral measure of A for a. 
As the integrand is bounded by 1 in modulus and tends pointwise to i, the 
dominated convergence theorem implies 



lim sh(si) = i\\a 



|2 



(Hi) =>- (iv): Obvious. 
(iv) => (i) By (8.3), 



2lmh(si) = 2is(v sh v si ). (8.10) 

By (8.8) and (8.10), there is a sequence s n such that — is n v Sn \ has a weak 
limit. Call this limit a. By (8.3) we have 



h(z)-h(s n \) = (z y v z ,v Sni ) + (v z ,-is n v Sni ). (8.11) 
Take the limit in (8.11) as s n — > oo to get (8.5). □ 

Theorem 8.12. Let h be in the Pick class of two variables, and assume 
h satisfies (8.8). There exists a representation as in (8.6) with 7-L finite 
dimensional if and only if h is rational and real-valued on the complement in 
M. 2 of its polar set. 

Proof. If h has a representation as in (8.6) with % <i-dimensional, it is clear 
that h is rational of degree at most d in each variable, and that h is real on 
M. 2 off its polar set. 
For the converse, let 

be a linear fractional map that maps the unit disk D to II, and 

0(z) = 

Z + I 

be its inverse. Let 

0(Ai,A 2 ) = 0oh(a(Xt),a(X 2 )). 
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This is a function in the unit ball of if°°(B 2 ), the space of bounded analytic 
functions on the bidisk. Moroever is rational if and only if h is, in which 
case they have the same bidegree, and is inner if and only if h is real- valued 



a.e. on M 2 . 



Assume h is rational and non-constant of bidegree (di,^)- By a result 
of G. Knese [10], there are Hilbert spaces Ai l and Ai 2 of dimension d\ and 
o?2 respectively, and analytic functions u 1 : D 2 — > Ai 1 and u 2 : ED 2 — > Ai 2 so 
that 



1-0(A)0(C) = (l-AxCiX^A),^)) + (1-A 2 C 2 )(« 2 (A), M 2 (0). (8.13) 
Define functions v r : II 2 — > Ai r for r = 1, 2 by 

Then an algebraic manipulation transforms (8.13) into 



h(z)-h(w) = (z 1 -w 1 )(v 1 (z),v 1 (w)} + (z 2 -W2){v 2 (z),v 2 (w)}. (8.14) 
Let 

K(z,w) = (v 1 (z),v 1 (w)) M i + (v 2 (z),v 2 (w)) M 2. 

This has rank less than or equal to d = di + d,2, so it is the reproducing kernel 
for some Hilbert function space Ai on II 2 of dimension less than or equal to 
d. By (iv) =^ (i) of Theorem 8.4, we have a vector a such that (8.5) holds. 
Now follow the proof of (i) =^ (ii), and observe that since B is defined on 
a finite dimensional space, its defect indices must match, and so it can be 
extended to a self-adjoint operator A on At. □ 
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